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Abstract 

The standard definition of tlie electromagnetic radius of a charged particle (in par- 
ticular the proton) is ambiguous once electromagnetic corrections are considered. We 
argue that a natural definition can be given within an effective field theory frame- 
work in terms of a matching coefficient. The definition of the neutron radius is also 
discussed. We elaborate on the effective field theory relevant for the hydrogen and 
muonic hydrogen, specially for the latter. We compute the hadronic corrections to the 
lamb shift (for the polarizability effects only with logarithmic accuracy) within heavy 
baryon effective theory. We find that they diverge in the inverse of the pion mass in 
the chiral limit. 
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1 Introduction 



The possibility to measure the muonic hydrogen (fip) Lamb shift at the 30 ppm level at the 
PSI [1] and, consequently, to obtain the proton radius with an improved accuracy of, at least, 
one order of magnitude better than current determinations, has triggered a lot of theoretical 
effort in order the theoretical errors to compete with the estimated experimental ones. We 
refer to [2, 3, 4] and references therein for a detailed explanation of the actual status of the 
art. On the other hand, it is also planed to study the muonic hydrogen properties at a future 
Neutrino Factory Complex [5]. 

Measuring the proton radius with such accuracy poses some theoretical questions on the 
definition of the proton radius. Actually, these problems already appear with the present 
claimed precision of around 1.5 % (see for instance [6]). The problem has to do with the fact 
that the proton radius is only defined at leading order in the electromagnetic coupling a. 
Beyond this order, the usual definition through the use of form factors runs into problems 
because it becomes infrared divergent. This fact has already been discussed by Pachucki 
[7]. Here we address the problem using potential NRQED (pNRQED) [8, 9, 10], an effective 
field theory suited for the study of non-relativistic systems. We will closely follow Ref. [11], 
where the hyperfine splitting of the muonic hydrogen and the hydrogen (ep) was studied (we 
refer to this reference for further details, definitions and notation). Within this framework, 
and beyond leading order in a, we will see that what one is really measuring is a matching 
coefficient of the effective theory, which is dependent on the scale (scheme) at which the 
matching has been performed. 

On the other hand we will see that there is a series of hadronic contributions that can be 
computed using heavy baryon effective theory (HBET) [12]. In particular, the non-analytic 
(and model- independent) dependence on the light quark masses of the Lamb shift can be 
obtained. We will do so here for the Zemach and vacuum polarization piece of the hadronic 
correction. The polarizability corrections will only be considered with logarithmic accuracy. 
The full computation will be presented elsewhere. 

2 Effective Field Theories 

We will obtain pNRQED by passing through different intermediate effective field theories 
after integrating out different degrees of freedom. The path that we will take is the following: 

HBET ^ (QED) ^ NRQED ^ pNRQED, 

much in the same way as we did in Ref. [11]. The fact that we are interested in the Lamb 
shift will force us to consider other operators which were not considered in that paper. 

2.1 HBET 

Our starting point is the SU(2) version of HBET coupled to leptons, where the delta is 
kept as an explicit degree of freedom. The degrees of freedom of this theory arc the proton, 
neutron and delta, for which the non-relativistic approximation can be taken, and pions, 
leptons (muons and electrons) and photons, which will be taken relativistic. This theory has 
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a cut-off n « nip, and much larger than any other scale in the problem. The Lagrangian 
can be split in several sectors. Most of them have already been extensively studied in the 
literature but some will be new. Moreover, the fact that some particles will only enter 
through loops, since only some specific final states are wanted, will simplify the problem. 
The Lagrangian can be written as an expansion in e and l/rrip. and can be structured as 
follows 

£^HBET — -\- Ci -\- Ct^ -\- Ci-^ + >^(JV,A) + ^{N,A)l + ^{N,A)Tr + ^{N,A)ln, (1) 

representing the different sectors of the theory. In particular, the A stands for the spin 3/2 
baryon multiplet (wc also use A = ttia — trip, the specific meaning in each case should be 
clear from the context). 

Our aim is to obtain the Lamb shift with 0{m^.a^ /m^F{miJmT^,mT^/ ^) accuracy. In 
particular, we want to obtain the non-analytic behavior in the masses of the light quarks 
and in the A. Actually, the full dependence on the ratio m7r/A can be obtained using chiral 
Lagrangians. In this paper we will mainly focus on a partial contribution to it: the so-called 
Zemach correction. 

To obtain this accuracy, we need, in principle, Eq. (1) with Oil/m^ accuracy. The 
different pieces of the Lagrangian can be found in Ref. [11]. There are a few differences due 
to the fact that we are now concerned with the Lamb shift correction. This implies that 
some terms in the nuclcon Lagrangian, which were not considered in Ref. [11] have to be 
considered here. In particular, the term 

where d2^R stands for the hadronic contribution. This term can be eliminated by field re- 
definitions (see for instance [13]) but we will keep it explicit since this contribution can be 
singled out by considering the lamb shift of a purely leptonic bound state. 

For the nucleon there is a new term that has to be considered (actually in this one is 
where the proton radius is encoded): 

5>CAr = A^;|-e^[V-E]|Ar^. (3) 

Specially relevant in our case is c||jj (see Eq. (14) in [11]). Nevertheless, the leading 
contribution in a (of 0{a^)) is suppressed by an extra factor mijuip (see Eq. (39)), i.e. 
c|}j ~ a^mijmp, which goes beyond the aimed accuracy of our calculation and we can 
neglect it. 



2.2 NRQED 

Since we want to obtain corrections to the lamb shift, besides the operators considered in 
Ref [11] new operators have to be considered with 0(l/(mpm7r)) precision. This means that 
some operators of up to dimension 7 have to be included in the Lagrangian. They read 
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16m^ 



>iV,. (4) 



c^j, correspond to the proton polarizabilites. The relations are the following: 



ct=Aml^, ct = -'^(af+f>t\ (5) 



a a 



We also note that the operators considered in Eqs. (2) and (3) also appear here. In 
the first case d2,R d2,NR- d2,NR is scale independent but pion effects should be included 
now. We will discuss in further detail these corrections in sec. 4. c^\i^) is a matching-scale 
dependent quantity and has already been considered in HBET. The difference is the scale v 
considered in each case, whereas in HBET ttIt^ <^ v <^ nip, we now have mi^a <^ v <^ m^. 

This means that in HBET, effects due to the pion scale were not included in c"i^{i') but now 
they are. A similar discussion applies to c^'jvii(^) (see Eq. (17) in Ref. [11]), which is also 
scale dependent. 

2.3 pNRQED 

After integrating out scales of 0{mi^a) one ends up in a Schrodinger-like formulation of the 
bound-state problem. We refer to [8, 9, 10] for details. The pNRQED Lagrangian for the ep 
(the non-equal mass case) can be found in Appendix B of [10] up to 0{ma^). The pNRQED 
Lagrangian for the np is similar except for the fact that light fermion (electron) effects have 
to be taken into accoimt. The explicit form of the Lagrangian has not been presented so far. 
We do so here. The explicit Lagrangian reads (up to 0{ma^)) 



LpNRQED = J d'^d'XS\^,X,t)l^ido - + ^ + 

-y (x, p, o-i, <T2) + ex • E(X, t) |>S(x, X, i) - y d'X^F^^F'^^ , 



8mf 2M 



(6) 



where M — -|- m^, //12 = m'^^ i ^ ^^<i X, and p and P are the relative and center of 
mass coordinate and momentum respectively. 

V can be written as an expansion in 1/m^, l/nip, a, ... We will assume 1/r ~ rrie (which 
is realistic for the case at hand) and concentrate in the l/m^^ expansion. We assume that 



<< nip. 



In order to reach the desired NNNLO accuracy in q;^, V^°^ has to be computed up to 0(q;^), 
V^^^ up to 0(0'^), y'-^-' up to O(tt^) and V^^^ up to 0{a). Nevertheless, there is not contribu- 
tion to V'^^^ at tree level. We first write the potential up to order 1/m^ in momentum space. 
Order 1/mP. We define V^*^°^ (the Fourier transform of V^^^): 

y(o) = -4TrZ^Zpav{k)^, (8) 



3 



which in fact defines ay, which is gauge invariant. There is another very popular definition 
that enjoys all the nice properties we would like a to have (gauge invariance,scheme/scale 
independence...): 

«^//(^) = °' i + n(-k2) ' 

where 

U(e) = aU^'\e) + a^U^^\e) + a'U^'\k'') + ... 

is the electron vacuum polarization and n(0) = 0. a^ff corresponds to Dyson summation. 
If we express av{k) in terms of aeff{k), we have 

av{k) ^ aeffik) + ^ Z^Z^^a^fik) ^ aeffik) + 5a{k) , 5a{k)^0{a*). (10) 

n,m=0 
n-\-7n=even>0 

The sum is constrained to fulfill n + m — even because of the Furry theorem. Each a^J^jP\k) 
is also gauge invariant. In order to achieve 0{ma^) accuracy we need to know H^^^ 

n(3) 

(see [14]) and the leading, non- vanishing, contributions to ^^^^■'(A;), o^^^ff^ and a^^j^K A 
discussion on the latter (which still remain to be computed) can be found in Ref. [15]. 

As we have already mentioned, the matching is performed in a expansion in 1/m, a and 
the energy. The leading-non-zero contribution to the potential from the expansion in the 
energy appears at 0{l/mP) and 0{a'^) and it will give rise to a contribution of 0{ma^). It 
follows from the electron vacuum polarization correction to the potential shown in Fig. 1. 
In momentum space it reads 

5VvP = Z^Zpe^^nik^) , (11) 

where 7r(/c^) = aU.'^^^k'^). We note that this correction depends on ko. Since ko ~ ma^ 
and k ~ ma we can expand in j^. The leading order is ko independent and produces the 
standard vacuum polarization correction that we have already taken into account in a^ff: 

Z,Z,e'^n{-k') . (12) 



The next correction reads 



6Ve = V§^'(-k') , (13) 



and is new. 

The explicit expression of tt reads 



TT J4 q^{mjq^ — k^) 



where 



^(0 = ^1-^(1 + ^ 



3 V q^ \ q^ 

The ko dependence in the potential translates into terms with time derivatives in the pN- 
RQED Lagrangian (or in other words an energy dependent potential). At the order we are 
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El, p 



ko, k 




Figure 1: Leading correction to the Coulomb potential due to the electron vacuum polarization 
in the Coulomb gauge, k = p — p' and ko = Ei — E[. 

working, one can get rid of these time derivatives by using the equations of motion. Note that 
we refer to the full equations of motion (including the potential), and not just the equations 
of motion of free particles, reflecting the fact that the matching has to be done off-shell. The 
procedure has been illustrated in the Appendix A of [10]. Effectively, this transforms SVe 
into the following expression 

We see that it becomes a 0(l/m^) potential. Field redefinitions should produce an 

equivalent result. We note that if wc want to go to higher orders than 0{ma^), we should 
work with the former expression Eq. (13), or work in a more general way by using field 
redefinitions. 

Order 1/m. They will give rise to 0{ma^) corrections so they will be neglected in the 
following. 

Order 1/m^. At the order we are working they can be written as follows^: 



naeffik) 



— 5" + ^11 — 5" 



(16) 



V-W = -i2.<,.„(*)^ ■ { K-^ + \ , (17) 




l>« = -Z,Z,16™(^^ + :^ + ^| , (18) 

4™^,^_u^, (19) 

rriijrnp \k.-^ k^ ) 

Vf'> = _':^EhlM(P^ . (z^eWs, + Z„c!?'s,) , (20) 



^Wc have slightly changed the definitions with respect [11]. We now set Zi, = Zp = 1 in the terms 
proportional to non-trivial matching coefficients in the Lagrangian (ie. Cf, cg, • • •). This allows to deduce 
the neutron case by setting ^ afterwards. 
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^^'^ = + 3<^i^) - 2<]vi.S^} . (22) 

nip 

The classification of the diagrams follows the one of Fig. 1 of Ref. [10] generahzed to the 
non-equal mass case and with the replacement a — > cte//- This ensures that the expressions 
for the potentials are correct at one- loop (the one- loop vacuum polarization is included). 
Written in this way subleading terms in a are also included. It should be keep in mind that 
at some point the replacement of a by a^ff is not enough and new corrections have to be 
included in the potentials. 

2 2 
Tn Tfl 

d2,NR = d,,n + -fKAO) = ^n;(0) , (23) 

where n'^(O) is the derivative of the hadronic vacuum polarization (we have defined n/j(— k^) = 
— k^n'^(O) -|- . . .) and n'^,r(0) the contribution due only to pions. 

The one loop diagrams of fig. 2 of Ref. [10] are equal in our case up to a trivial mass 
rescaling: 

The procedure we have followed to obtain the potentials is different from the one used 
in Ref. [2]. For instance, in this last reference no mention to a correction to the type of 
Eq. (15) is made. To facilitate the comparison one should work in position space (it is also 
eventually useful to write the potential in a basis where the angular momentum structure of 
the potential is more visible) . The non-trivial comparison can be reduced to the contribution 
due to SVe and the one- loop correction to V^^\ which in our case reads 

Hoop TT mfjUp J4 [meq) I 2 I r V 2 



'^^'^ 1 + meqr) +i!M)!e-e.. (^^ + !ML^ _ 2^S{r) \ , (26) 



2r3 



where L is the angular momentum. We have checked that this expression agrees with the 
corresponding expression in Ref. [2]. 

In this paper, our main focus are the hadronic corrections, in particular those propor- 
tional to the 5{r) potential^. They can be encoded in the matching coefficient of the spin- 

^There are other hadronic corrections. In particular, one also has to consider the one-loop vacuum- 
polarization correction in Eq. (16), proportional to c^-^\ Nevertheless, this term has a different functionality 
on r, which makes it possible (at least as a gedanken experiment) to single it out by considering different 
observables (the result will depend on m^a/(nme) for the spectrum or on |p|/me above threshold). Therefore, 
we will not consider it in this paper although it should be considered for the Lamb shift of low excitations 
of muonic hydrogen with 0{ma^) accuracy (see [2] for an explicit computation). 
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independent delta potential: 



7~)had 

5V = :li^^A(3) 



mt 



where 



(27) 



(28) 



The spin-dependent term does not contribute to the average energy over polarizations. 
Therefore the hadronic contribution to the lamb shift reads 



5E = {E{s) - E{p)) 



7~)had 

ml 



1 (ii^pa 

OlQ— I 



(29) 



3 Definition of the proton radius 

The usual definition of the electromagnetic radius of the proton is the following 



(30) 



q2=0 



This definition runs into problem once electromagnetic corrections are included in the com- 
putation of Gp^E because it becomes infrared divergent (this is so for a charged particle, for 
the neutron this definition is infrared safe). This problem has already been pointed out in 
Ref. [2]. Nevertheless, this is not a problem from the effective theory point of view, where 
the parameters of the theory (the matching coefficients) do depend in general on the scheme 
and scale on which the computation has been performed. Therefore, our definition of the 
proton radius is the following 

» - 4 



1 - 2 2 

1 = 3^p"^p > 



which naturally relates to the one in Eq. (30), since 



cd = 1 + 2F2 + 8F[ = 1 



5m, 



dq^ 



(31) 



(32) 



q2=0 



where we have already used that Fi = 1 due to gauge invariance. Eq. (31) is the natural 
definition from the point of view of effective field theory (a matching coefficient). Another 
issue is the separation of "hadronic" from electromagnetic corrections. We stress that the 
electromagnetic effects are included in the form factors above. We note that Tp is now an 
scale dependent quantity and fulfils an expansion in a: 



(33) 



We will define the pure hadronic definition of the radius ("hadronic proton electromagnetic 
radius") to be the first term in this expansion, Note that this quantity is not an 

observable from the effective field theory point of view, since it will always appear (in the 



7 



same combination) with another hadronic corrections in the observable. This definition for 

the "hadronic" piece is different in general^ of the definition in [16]. 

We can also know what the running of each term of Vp is. As we have already mentioned, 
the first term is scale independent but not the second (even if this running can be understood 
in terms of point-like particles) . 

^;frf = 0, r^vA,(i) = _li^, .... (34) 

In summary, if we give a number for the proton radius we have to say in which scheme (for 
instance MS) and scale it has been obtained if a precision of around 1 % is claimed. In Ref. 
[7] a prescription has been taken to fix the ambiguity with 0{a) precision (how to generalize 
this prescription to higher orders in a is left unsettled). There, the scale dependence at 0{a) 
is regulated by the proton mass and the finite pieces fixed by performing the computation 
assuming the proton to be point-like at scales of order nip. Nevertheless, this assumption 
is not justified since at these scales several resonances appear. We would also like to stress 
that our definition of the proton radius through a matching coefficient, c^£\ makes easier to 
obtain the dependence of any observable in this quantity, in particular for the cross section 
a{e^p e~p). 

What we have said is quite general for any nucleus. Indeed, the object which is really 
measured at very low momentum transfer (above or below threshold) between the nucleus 
and lepton is Dd{u) = + D]^^'^. Therefore, one could think that it is more proper to 
give numbers for this object (for a given scale and scheme). Nevertheless, if we want to 
get some understanding on the structure of the nucleus it maybe convenient to distinguish 
the different hadronic and leptonic contributions of D^iy) as some of them can be singled 
out from different measurements. For instance the contribution due to d2,NR, as well as 
the leptonic contributions, can be singled out since they can be computed explicitly and 
measured changing the nucleus by another lepton. Moreover, it also makes some sense to 
distinguish between c|'ivR ^^'^ '^d (^^ least at 0{a'^mi.)), since the leading dependence on 
c^']vR is proportional to the lepton mass and can be distinguished by changing the lepton. 
Beyond or beyond m;., this distinction is not that clear and can not be achieved in general. 
This has to do with the fact that we are not working in a minimal basis of operators in the 
effective Lagrangian. In other words, one can always perform a field redefinition to eliminate 
c'^\ which gets absorbed in c^^nr- Therefore, in the long term one should work with D^'^'^, 
which can be related with observable quantities (for the neutron case it can easily be related 
with the scattering lengths"^). 



■^To get a closer connection between both definitions one should take /Ui < nie in Ref. [16], where /xi 
is the matching scale at which the parameters of the theory are make equal between the hadronic and 
hadronic+clcctromagnctic case. In other words, one natural place to make the connection between hadronic 
and hadronic+electromagnetic quantities is at very low energies (/Ui 0), where QED has an infrared fixed 
point and consequently the fundamental parameters of QCD do not run. 

^The definition of the neutron radius could also be worked out on similar terms. One could write the 
effective Lagrangian, which would be similar to the proton one but for a particle with zero charge, i.e. with 
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4 Matching HBET to NRQED 

The hadronic corrections to the spin-independent delta potential are encoded in the correc- 
tions to d2,NR, and c^^nr- 

We will only compute the hadronic corrections to d2,NR and cz,nr (due to energies of 
0(m^, A)), since our aim is to obtain c^^ from experiment (for instance, from the muonic 
hydrogen lamb shift). If we were interested in r^^ one should also compute the hadronic 

0{a) corrections to c£\ 

The leading correction to d2^NR is due to the one-loop pion correction to the vacuum 
polarization. The diagram that contributes at leading order is the same than in Fig. 1 
changing the electron by a pion in the limit of transfer momentum going to zero. It reads 

The contribution to c^,^nr (sec Fig. 2) can be written in a compact way in terms of 
the form factors of the forward virtual-photon Compton tensor (see [11] for details). It 
reads (this expression is only valid in 4 dimensions and in D-dimensions has to be properly 
generalized) 

c^^NR = i9'm,m,^ j (^I^F^^pf {^^^^°' + " ^'^'^^'^ ^'^} " ^^^^ 

This result keeps the complete dependence on m/.. This contribution is usually organized in 
the following way 



the replacement of covariant by normal derivatives. The definition of the neutron radius is also 

(35) 



r: 



2 a dGn,E{q^) 



and from the effective field theory point of view 

4")^^« = 2Ff)+8F(")', (36) 

where we have already used that F^"^ = due to gauge invariance. The point that we would like to remark 
here is that rj; does not correspond to the neutron scattering length once 0{a) corrections are included. 
Rather the neutron scattering length 6„;. (where li represents the specific lepton) reads (note that d2,NR 
does not contribute in this case) 



bnu = -r— ^Zx,ct^ - -h:nb. : (37) 



1 7 (n) 2 

This is a low-energy constant (formally we can take the limit of transfer momentum going to zero) and 
measurable quantity which can (in principle) be obtained with arbitrary accuracy from experiment. Note 
also that if one works in a minimal basis of operators in the effective theory (by absorbing in Cg'j^^j 
through field redefinitions) there is a one-to-one correspondence of the scattering length with a matching 
coefficient of the effective theory: Cg'^^. 
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Figure 2: Symbolic representation of Eq. (39). 

We remind that is O^a^mijmp) and can be neglected. The second term corresponds to 
the diagram assuming the proton non-relativistic and point-hke. With the precision needed 
it reads (in the MS scheme) 




Figure 3: Symbolic representation (plus permutations) of the Zemach correction in Eq. (42). 

The second term in Eq. (40) is the so-called Zemach correction, which is one of the main 
concerns of this paper. It is symbolically pictured in Fig. 3 and reads 

S4%emacH = 4(4vra)2mjm, / ^^.^.G^e^G^ . (42) 

This term appears to be finite to the order of interest (it produces no logs) and it agrees 
with the result obtained by Pachucki [2] at leading order. Note, that this result holds for 
both ep and /ip. In other words, the exact dependence on mi- is kept. gP = 1. We take the 
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expression for from [17]. The use of effective field tfieories and dimensional regularization 
is a strong simplification: we only need the non-analytic behavior of G^^ in g^, the analytical 

behavior produces scalelcss integrals, which are zero in dimensional regularization. This is a 
reflection of the factorization of the different scales. We do not need to introduce the point- 
like interactions to regulate the infrared divergences of the integrals at zero momentum. A 
similar thing happened [11] for the computation of the hyperfine splitting of the hydrogen, 
where, again, there was not not need to introduce the point-like interactions to regulate the 
infrared divergences of the integrals at zero momentum. 
The result of the computation reads 




where 

B^+2r - -f^B,+2r > , T > , (44) 

and the definition of Bn and Hn can be found in Appendix A. The term proportional to 
are due to scales of 0(A), the contributions proportional to Hn are due to scales of 0(m7r). 
We see that the latter give contributions to the Lamb shift which are non-analytical in the 
quark mass. 

We have then obtained the leading correction to the Zemach term in the chiral count- 
ing (supplemented with a large Nc counting). This is a model independent result. Other 
contributions to the Zemach term are suppressed in the chiral counting. 



(-l)'T(-3/2) f 
r(r + l)r(-3/2-r) \ 




Figure 4: Diagram contributing to the polarizability correction with Inmg accuracy. The 
matching coefficients of the proton can be ca^ or ca^, or, in other words, the proton polariz- 
abilities. 

Finally, we consider the polarizability corrections. For those we would need the full 

expression for the virtual-photon Compton tensor at one-loop in HBET, which is lacking at 
present. This computation is relegated to a future publication. We will only consider here 
the limit m^. -C m7r, which is relevant for the hydrogen atom, and only consider the logarithm 
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(In me) enhanced terms. These can be obtained by computing the ultraviolet behavior of 
the diagram in Fig. 4. This contribution is proportional to c^j and 0^2- It can be written 
in term of the polarizabilites of the proton (see [18, 19]) 



Pm 



ln{mi.) . 



(45) 



For the pure pion cloud, the polarizabilities were computed in Ref. [20]. The contribution 
due to the A can be found in Ref. [21]. The scale in the logarithm is compensated by the 
next scale of the problem, which can be m.,^ or A. For contributions which are only due to the 
A or pions, the scale is unambigous. In the case where pions and A are both present in the 
loop we will choose the pion mass (the difference being beyond the logarithmic accuracy). 
The result reads 



2 2"^ii.2 1 A 49 2 2 
^W^,ln — + -na 9a- 
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mi 



In(^) 
mi 



mt 



45A 



mi 



(46) 



AA^ -Adml , „ \ , ,m^, 

+ ^2 ^Ini? In — , 

A^ — mi I mi 



where 



mi 



(47) 



We would like to stress here that some of the corrections considered in Eq. (46) are missing 
in the dipole approximation [22], which is sometimes used in the literature. For instance, the 
first term in Eq. (46) (the pure Delta correction) would be missing in that approximation. 



5 Phenomenological analysis 

Our aim here is to see how much of the hadronic correction can be understood in terms of 

the physics of the m^^ and/or A scales. 

There have been some previous estimates [2, 3, 23] of the Zemach correction using phe- 
nomenological model parameterizations of the proton form factors. They obtain quite similar 
numbers. For instance, using the dipole parameterization of the proton form factor, 

G^(q^) = , ^ 2,2 , (48) 

H) 

Pachucki [2] obtains 0.021 meV (slightly larger in [3] or if one takes into account relativistic 
effects) using A = 848.5 MeV. Within chiral perturbation theory, the Zemach correction to 
the muonic hydrogen Lamb shift in the limit A — > cxo reads 

= -0.01008 meV. (49) 



^^Zemach^^ = 2) 
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This is around 1/2 the standard predictions using the phenomenological parameterization of 
the form factors. The inclusion of the effects due to the A increase this number by around 
a factor of two and bring it surprisingly close to the phenomenological estimates. 



^^zemach^^ = 2) = -0.01922 meV . (50) 

The corrections due to the Delta are large even if the series in m-j^j IS. is convergent (—0.01008— 

0.00561 H ). It would be interesting to study whether something similar happens for the 

hyperfine splitting. 

Obviously, for the hydrogen we also obtain similar numbers to the phenomenological 
estimates if the A is included 

^^Zemach^^ = 2) = -3.49652 Hz . (51) 



The pure pionic correction reads 

^^Zemach/ _ 3) = -1.83404 Hz . (52) 

^ ' A— >oo ^ ' 

The polarizability correction shown in sec. 4 is only reliable for the case of hydrogen. It 
is proportional to the electric and magnetic proton polarizabilities. Therefore, the question 
is how good HBET can reproduce the experimental number for the polarizabilities. Actually, 
if only the pure "pionic" contribution is included a very good agreement is obtained. This 
agreement is deteriorated if corrections due to the A are included (specially for /9m)- In any 
case, the correction we obtain for the Lamb shift correction reads 

64.4841 



V? 



Hz , (53) 



and including the A 

^-r^-^Hz. (54) 

We see that they are reasonable close to each other and with the value obtained using the 
experimental number for the proton polarizabilites [19]. Slightly larger numbers are obtained 
using the experimental data of the deep inelastic structure functions [23] or models [24] . 

Finally, wc would like to stress that in our determination there is no free parameter. The 
leading contribution is the one we have computed. This is due to the chiral enhancement. 
Our result diverges in the chiral and/or 1/iVc limit. 

For the vacuum polarization, the contribution due to the pion loop is numerically small 
(even if leading from the power counting point of view). One obtains njj ,^ ~ l.Ox 10^"^ GeV^^. 
This should be compared with the experimental figure for the total hadronic contribution 
H^ ~ 9.3 X 10~^GeV~^ [25]. Actually, the main contribution comes from the p even if it 
is suppressed in the power counting (~ 6.1 X lO-^GeV"^). In any case we need H'^^ with 
(a / rrip) precision for which wc can take the experimental number for the vacuum polarization 
(matching coefficient di^NR) and the correction to the Lamb shift reads 

5^^^. pel. ^ JJ^^ kHz , bE-^- = -^^^ meV , (55) 
which for the hydrogen case reduces to the result obtained in Ref. 
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6 Conclusions 



We have studied the lamb shift of hydrogen and muonic hydrogen within an effective field 
theory framework, connecting the physics at the hadronic scale (HBET) with the physics 
at the atomic scale (pNRQED). We have focused on the hadronic corrections and showed 
that a natural definition of the proton radius can be given as a matching coefficient of the 
effective theory. 

The leading hadronic corrections to the Lamb shift are a prediction of HBET (unknown 
couterterms are subleading). We have exactly computed them for the Zemach and vac- 
uum polarizability terms, and for the polarizability term with logarithmic accuracy. For 
the Hydrogen atom the computed correction is beyond the present experimental accuracy. 
Nevertheless, the result is perfectly relevant for the muonic hydrogen, for which it would be 
possible to get a model independent determination of the proton radius if the polarizability 
corrections were computed in HBET. 

We have also done some partial checks of the (QED) hydrogen muonic computation. 
Nevertheless, it would be useful to perform a comprehensive numerical study for muonic 
hydrogen (since some pieces of the QED computation has only been performed by one 
group). Moreover, some further work remains to be done for the hadronic corrections: 1) 
to obtain the complete result for the polarizability correction to the muonic hydrogen Lamb 
shift and 2) to perform the same analysis for SU(3) (including kaons). 
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Appendix A 



I" 

Bn = / 
^0 



dt , = In 



(56) 



v^r(; 



n 



n. 



I HarmonicNumber( ) — HarmonicNumber(l 



4r(2-f) V 

, oi-n, r(n-2) 1 n-2 n-l _n n 

+ 2 TT 3-^21-, ^ , ^ io^oi^) 



r2(|) ^ ^4' 2 ' 2 ' 2' 2^ 



-2+n 



TX 



r(2-f) 



-3 (n-2) r(l -n) + 2r(-n) 



r(3 - n) 

{n — 3) n 
2 d 



( 3 n 

-3 + (n — 3) n ( tt cot(n7r) — log(2) — ■^(- — — ) + '>\){n) 



n — 3 dx 



2^i(i, — 7, — ; + ^; -1) 



a;=0 



where '^{n) — dr{n)/dn and r(n) is the Euler F function. 

^ n!(2n- l)!!r[-3/2] 
2(2n)!!F[l/2 + n] 



(57) 
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